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1 Introduction 

During inflation, all light scalar fields are subject to stochastic fluctuations [1-3]. Their 
amplitudes tend to increase by a process that could be described as random walk. Once 
inflation is over, the fields can start moving classically, with many interesting cosmological 
consequences. An example is Affleck-Dine baryogenesis [4], where the motion of a complex 
scalar field, initially displaced from the minimum of its potential by inflationary fluctuations, 
dynamically generates a baryon asymmetry. However, if there are several scalars that are 
coupled to each other, their dynamics during inflation can be rather complicated as the 
fluctuation of one scalar can backreact on the fluctuations of another. This is also the 
situation one would encounter in theories with flat directions (or moduli fields), which are 
rays in field space along which the potential either vanishes completely or is very small. An 
example is the Minimal Supersymmetric Standard Model (MSSM), which has some three 
hundred flat directions; these have all been classified (for a review, see [5]). In the MSSM 
flat directions correspond to configurations where some of the field values are related to each 
other while the rest are set to zero. As a consequence, most of the flat directions are not 
simultaneously flat, and if there is a large field fluctuation along some given flat direction, 
many other potentially flat directions will no longer remain flat. In effect, fluctuations may 
provide some of the flat directions with effective masses that prevent the growth of the 
amplitudes, as was discussed in [6]. 

The dynamics of such a system of coupled, stochastically fluctuating scalar fields is in 
general very complicated. In the present paper we attempt to give a careful analysis of the 
growth of fluctuations in a system that only has two scalars, coupled to each other. Unlike in 
[6] , we do not require that the initial field (expectation) values are zero but allow also for the 
possibility that at the onset of inflation, the light field (flat direction) can have an arbitrary 
expectation value. We assume that inflation is independent of these scalars, which for the 



- 1 - 



purposes of inflationary dynamics are just spectators, and that the Hubble rate H can be 
taken as a constant. We further assume an FRW universe with H = a/ a, where a(t) is the 
scale factor, and the number of e-foldings given is by N = ln(a(i)/a(0)). By definition, for a 
given background value <fio, the field 4> is called light if V"(0o) ^ H 2 . 

The paper is organized as follows. In section 2, the basics of the stochastic approach to 
inflation are reviewed. In section 3, we study the behavior of a massless scalar field coupled 
to a massive one, and approximate the behavior of the massless field when the former begins 
with a large expectation value. In section 4 we investigate the consequences of such an 
assumption and derive an expression for the inflationary fluctuations of the coupled massless 
field. In section 5 we present our conclusions. 



2 Stochastic formalism 



In general, the equation of motion for a minimally coupled, self-interacting scalar field with 
mass M and interaction potential V(0) is 

$ + + M 2 $ 2 - a~ 2 V 2 $ + V'($) = . (2.1) 

Provided M 2 <C H 2 , the field $ undergoes quantum particle production from its vacuum 
fluctuation during inflation. Owing to the probabilistic nature of the vacuum fluctuation, 
the process of particle production is stochastic. Also, due to the inflationary expansion, 
the momenta of the created quanta are redshifted as the corresponding field modes become 
superhorizon. Consequently, the evolution of the light field $ on long wavelength scales 
contains an stochastic component due to the inflationary particle production. 

In order to study the stochastic evolution of a light field <$ during inflation, one intro- 
duces a coarse-graining scale k s (t) = ea(t)H, where e <C 1 [1-3]. This scale serves as a divide 
to decompose the field $ into its long and short wavelength parts as follows 

$(t,x) = <f>(t,x) + <l> q (t,x), (2.2) 

where the short wavelength part of the field is given by 

/rl 3 k r i i 

^3 9{k - k s ) [a k <p k (t) e ik x + a{^ k (t)e- ik - x ] , (2.3) 

where k = \k\, describes the quantum behavior of the scalar field on subhorizon scales. 
The long wavelength part <j> (left unspecified), describes the field's classical behavior on 
superhorizon scales. 

The equation of motion of a scalar field in de Sitter space can be written as [7] 

+ 3iI0-a- 2 V 2 + y / (0) =f(s,t), (2.4) 

where 

£(x, t) = - (4> q + 3Hj> q - a- 2 V 2 4> q + V'{4> q )) (2.5) 

is a stochastic source of white source, i.e. (£(x,t)) = 0, where () stands for the vacuum 
expectation value of quantum operators. 

To simplify the discussion, we adopt the so-called separate Universe approach [8], ac- 
cording to which the Universe can be regarded as an ensemble of patches of size k' 1 in which 
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the field is approximately homogeneous and evolves independently of the other patches. Un- 
der this assumption, the term V 2 can be neglected in Eq. (2.4) and the stochastic source 
£ becomes space-independent within each coarse-graining volume. In the case of minimally 
coupled, light scalar fields, the self-correlation function of the stochastic source is 

<«*)f(0> = ^(t-O- (2-6) 

Regarding the above average ( ) , we wish to emphasize that for a free scalar field with mass 
m <C H, the quantum average (0 2 ) computed by field theoretical methods [9] agrees with 
the average over an ensemble of independent universes [10]. Therefore, when applied to a 
light scalar field tfi during inflation, we are entitled to interpret the quantum average ((j) 2 } as 
an ensemble average over coarse- graining patches. 

We want to address the growth of field fluctuations in the regime in which (f> performs 
a slow-roll motion due to the flatness of its potential. In such case, the solution to Eq. (2.4) 
can be approximated by neglecting the second time derivative (p. The behavior of the coarse- 
grained field is then determined by the Langevin equation 

+ V'((f>) ~ . (2.7) 

From the standpoint of stochastic equations, it is straightforward to write down the Fokker- 
Planck equation [11] 

dtPtf, t) = P(d>, t)) + ^ t) (2.8) 



3H ^' 7 8vr 2 

describing the evolution of the probability density that the field $ has an expectation value 
on a separate Universe, or coarse-graining patch. The first term on the righthand side 
accounts for the deterministic motion of the field, or drift term, and the second for its 
random motion. 

In the simple case of a free scalar field with mass m and the potential V((j>) = \ m 2 (f> 2 , 
whenever the field is displaced from the origin, the restoring force driving the field to <j) = 
increases proportionally to (j). For sufficiently large <j>, the restoring force compensates for the 
diffusive motion of the field and an equilibrium situation arises. Since in the equilibrium the 
probability distribution becomes stationary with dfP((J),t) = 0, demanding that P(<p,t) van- 
ishes sufficiently fast at cp — > ±oo one readily obtains the well-known equilibrium amplitude 
of the field fluctuations 

<*%=g£- < 2 - 9 > 

In relation to this result, in Ref. [6] it is shown that flat directions of the MSSM reach 
a stationary fluctuation amplitude similar to the above when they begin with a vanishing 
expectation value. This results arises because the coupling of the flat direction to other 
fields provides the former with an effective mass, which then results in the emergence of an 
stationary distribution. As we shall show in the following sections, if the coupled massless 
field (or the would-be flat direction) begins with a large expectation value (</>) 3> H, the 
amplitude of the field fluctuations reaches the equilibrium value only after a sufficiently 
long wait. However, the stochastic dynamics in the coupled massless field in this case differs 
substantially from the one studied in [6] . Also, and in agreement with recent numerical results 
[12], we find that the coupled massless field is able to fluctuate as if free. Nevertheless, in 



- 3 - 



our approach we discover that this regime of quasi-free fluctuations applies for a limited 
time only. Beyond certain timescale, the system engages into a lengthy process ending up 
with the massless field fluctuating with its equilibrium amplitude in the entire ensemble of 
coarse-graining patches. 

3 Behavior in the large field limit 

In this section we investigate a system of two scalar fields, denoted by <f> and x, minimally 
coupled to gravity during inflation. We assume that the energy density of both fields re- 
mains always subdominant, hence their evolution does not affect the inflationary expansion. 
Assuming that <\> is massless while x has a squared mass m 2 , the Lagrangian of the system is 

C = ±d„$d»$ + l -d^X ~ \rn\x 2 ~ V($, x) , (3.1) 
where the interaction potential is 

^>X) = ^Y, (3.2) 

and g is a coupling constant. Such an interaction potential is ubiquitous in quantum field 
theory, and its consequences have been extensively studied in the theory of reheating and 
preheating [13-16]. Moreover, this coupling has been shown to provide a mechanism to trap 
string moduli fields at points of enhanced symmetry [17, 18], also to result in trapped inflation 
[19-21], and to give rise to a source of non-gaussianity in the inflaton's perturbation spectrum 
[22-24]. Typically, the moduli field playing the role of the inflaton approaches a location in 
field space where the x field becomes massless. If the light field <1> becomes sufficiently close 
to the point of enhanced symmetry, a x-P ar ticle production mechanism is triggered which 
modifies the dynamics of <I>. In this paper we do not study this situation in any significant 
detail, but will invoke the existence and implications of such a process. 
The equation of motion for the fields $ and x is 

<E> + 3H<t> - a" 2 V 2 $ + g 2 x 2 ^ = , (3.3) 

X + 3H X -a- 2 V 2 x+(m 2 x + g 2 $ 2 )x = 0. (3.4) 
The effective masses of $ and x are determined by 

2 _ 2 2 2 - 2 , 2,t.2 /o r\ 

m^ = g x , m x = m x + 9 * ■ ( 3 - 5 ) 

To solve the dynamics of the long wavelength part of <£, denoted by <j> (see Eq. (2.2)), our 
starting assumption is that within coarse-graining patches $ is approximately homogeneous 
(hence determined by <p) and that it has an initial non- vanishing expectation value <po, i.e. 

$(x,0)«<£o, (3.6) 

large enough so that the field x is initially heavy, namely 

m x (0) 2 ~ 5 2 ^»F 2 . (3.7) 

Later on we consider and discuss the case for which rh 2 ~ g 2 & 2 - 
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3.1 Integrated coarse-grained dynamics 

Since % is a heavy field, it performs oscillations around x = in a timescale that is very 
short compared to the Hubble time. On the other hand, we wish to consider the situation 
in which the massless field $ receives a small effective mass due to its coupling to x> and 
hence performs a slow-roll motion. In such case, evolves in a timescale large compared 
to the Hubble time. This timescale separation between the evolutions of <1> and x implies 
that we can integrate out the field x an d compute the resulting effective equation of motion 
for the long wavelength part (p. Neglecting dissipation effects, which we justify later on, the 
equation for the coarse-grained light field <j) can be written as 

3H(j> + m\4> ~ f (t) , (3.8) 

where m 2 = {m\) = g 2 (x 2 ) and we have neglected <\> and the gradient a _2 V 2 (/>. Of course, to 
use (3.8) consistently we must secure that 4> behaves as a light field, namely that g 2 (x 2 ) *C H 2 . 
Negligible dissipation effects also imply that we can treat £(i) as a stochastic source of 
gaussian noise whose correlation function is given by Eq. (2.6). 

To compute the variance (x 2 ) on the coarse-graining scale k s we need to solve for the 
mode equation for the Xk modes 

Xk + 3H Xk + (™ 2 + ^ ) Xk = • (3.9) 

Owing to the slow-roll motion of <p we may assume that m 2 ~ g 2 4> 2 does not change sub- 
stantially from the time when k 2 /a 2 = m 2 until the fe-mode exits the coarse-graining scale. 
In such case, the solution to the mode equation (3.9) matching the solution Xk = -^7= e ik l aH 



in the limit k/aH — > 00 is approximated by 

Xk® = a^ 2 ^e^m H ^{k/aH) , (3.10) 

where v 2 = Q/A — m^/H 2 and is the Hankel function of first kind. Expanding the Hankel 
function for k <C aH we obtain 

\Xk\ 2 = ^M l) (k/aH)\ 2 



4a? H 

e -WW n [1 + coth(|^|vr)] 2 sinh(|^|vr) _ 1 + coth (\u\ir) 



Aa 3 H \u\tt Aa?H\i 



(3.11) 



In order to compute (x 2 ) on the coarse-graining scale we follow the regularization scheme 
by Vilenkin [25] (see also [26]). The computation of the quantum average (x 2 ) in such case 
boils down to integrating out the Xk modes with momentum below the cutoff k = He Ht . 
Since we have introduced the coarse-graining scale k s for the massless field cf>, we need to 
compute (x 2 ) on that scale in order to substitute in Eq. (3.8). Therefore, we multiply by 
6(k s — k) in the momentum integral to filter the appropriate superhorizon modes. Further 
neglecting the contribution to (x 2 ) from the pre-inflationary epoch ((x 2 (0))) we have 

(X 2 )IR=/ TTT,0(k s -k)\xk(t)\ 2 . (3.12) 



k=H 
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Using Eq. (3.11) we obtain 

where = ff _l loge _l is the time necessary for horizon-size modes to exceed the coarse- 
graining scale. Since t& is typically a few Hubble times, in the following we consider that 
t S> td at all relevant times. Also, using \u\ ~ m x /H, one finds that Eq. (3. 13) simplifies to 

(x 2 W)iR-^(^) 2 — £ 3 - (3.14) 

An equivalent manner to get to this result is by using the perturbation spectrum of a heavy 
field [27] 

to compute (x 2 )ir = / P\(k) ^- Owing to the scale-dependence, this equivalent approach 
leaves manifest that the main contribution to (x 2 )ir stems from modes with wavelengths 
comparable to the coarse-graining scale, whereas the contribution from modes with larger 
wavelengths is suppressed. It is then clear that the scale-dependence featured in Eq. (3.15) 
is the reason of the factor e 3 in Eq. (3.14). 

From Eq. (3.14) it follows that to 2 . = g 2 (x 2 )lR ^ H 2 , thus allowing us to neglect (j> in 
Eq. (3.8). Also, in the limit m x <C g4> we have (x 2 )ir °c hence quantum fluctuations of 

the x field give rise to a slightly tilted linear effective potential for <fi. Quantum fluctuations 
of the $ field also couple to Xi which could result in dissipation effects further altering the 
dynamics of <f>. We discuss this possibility in the next section. 

3.2 On dissipation effects 

To discuss dissipation effects we bear in mind that the fields $ and x are to be eventually 
identified with fields in a supersymmetric theory, in the context of which dissipation effects 
have been extensively studied (see Ref. [28] and references therein). In this section we thus 
regard the system in Eq. (3.1) as contained in the scalar sector of a more general supersym- 
metric Lagrangian supporting the dissipation mechanisms studied in the literature. In this 
sense, a dissipation two-staged mechanism commonly considered in the literature couples 
the inflaton field to a intermediate heavy scalar field x- I n turn, the heavy scalar x decays 
into light degrees of freedom at a rate larger than the expansion rate. To give an example, 
in the context of supersymmetry such an interaction structure can be accounted for by the 
superpotential [28] 

W = g$X 2 + hXY 2 , (3.16) 

where X and Y are superfields and (f>, x> an d y refer to their bosonic components. During 
inflation, the field y and its fermionic partner y remain massless, whereas the field x an d its 
fermion partner ip x obtain their masses through their couplings to (ft, namely m^ x = m x = gcf). 
In the low-temperature regime, defined by m x > T, the dissipation coefficient is of the form 

T~CVT 3 /</> 2 , (3.17) 

where ~ h A N x Nl ccay , M x is the multiplicity of the X superfield and A/"decay is the number 
of decay channels available in X's decay. Although in our case the field <j> does not play the 
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role of the inflaton, the dissipation mechanism just described equally applies to (f) since in 
the regime of interest to us x is a heavy field. 

Dissipation effects alter the dynamics of the coarse-grained field in the strong dissipation 
regime, defined by the condition T > 3H. Using the slow-roll equations of motion for (j> and 
the radiation density [28] and the effective potential V e s(4>) = ^ <? 2 (x 2 )iR0 2 > it can be shown 
that the condition T > 3H translates into the bound 3> (ge 3 )~ 3 / 2 (cf)/H) 2 . Taking typical 
values of g and e and using <j) > H, the resulting value for is very large indeed. We thus 
conclude that dissipation effects can be safely neglected unless the x held belongs to a large 
representation of some Grand Unified Theory model [29]. Owing to the coupling between 
4> and the radiation bath, dissipation effects can also affect the perturbation spectrum of (f> 
even in the weak regime, defined by T < 3H. However, using again the slow-roll equations 
of motion for the field/radiation system it follows that the condition p r > H A (necessary 
for dissipative effects to alter the perturbation spectrum of (j>) translates into a bound on 
Cfi much tighter than the former one. In turn, this demands that the x held belongs to 
a representation much larger than those considered in Ref. [29]. Therefore, we disregard 
dissipation effects on the perturbation spectrum of 4> as well, which then allows us to use the 
correlation function in Eq. (2.6). 

It is worth clarifying that dissipation effects are so suppressed in this case because their 
magnitude depends on the kinetic density of <f>, which in turn depends on the steepness of 
T4fr (</>). Consequently, given the extreme flatness of the effective potential, dissipation effects 
can only become important when the smallness of the kinetic density of cj> is counteracted by 
a large number of degrees of freedom taking part in the dissipation mechanism. 

3.3 Effective Fokker-Planck equation 

After computing (x 2 )m we are ready to discuss the stochastic dynamics of long wavelength 
field (p. Although Eq. (3.8) can be readily integrated to work out the two-point correlator 
(cj) 2 ), we find it more convenient for our purposes to write down and solve the Fokker-Planck 
equation (2.8). Using V^{4>) = g 2 (x 2 )iR4> with (x 2 )ir as given by Eq. (3.14) we have 

g 2 H 2 e 3 n ( <f>P((f>,t) \ H 3 « , 
dtP = y ~^r d* -f - )+^d 2 P. (3.18) 



36vr 2 v yifnl+g 2 ^ 2 ) 1 / 2 

We wish to find an approximate solution to this equation for m 2 <C g 2 4> 2 - Restricting our- 
selves to 4> > 0, the force provided by the effective potential is determined by 

VM) gH 2 e 3 , s 

if^W' (3 - 19) 

which is constant. The relevant question to address is whether such a constant force driving 
the field to = manages to restrain the field's diffusive motion. The answer to this question 
is straightforward: since the effective potential is a confining one in the entire domain of field 
values, the probability distribution becomes stationary after a sufficiently long wait [30]. 
Although the approach to equilibrium in a linear potential can be studied analytically using 
the method of eigenfunctions [30], a simpler approach to the problem is possible because in 
the regime of interest to us (0 < fh x <C g4>) the Fokker-Planck equation can be integrated 
directly (equilibration timescales for the case of single field has recently been investigated in 
this context in [31]). Neglecting the term in m 2 in Eq. (3.18) we have 
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We emphasize that since this equation holds for < m x <C g4> only, we should not expect to 
recover the equilibrium distribution that follows from imposing dfP = in Eq. (3.18) (valid 
for all 0) by taking the limit t — > oo in the solution to Eq. (3.20). Nevertheless, as explained 
in the next section, the solution computed below is sufficient for our purposes. 

To integrate Eq. (3.20) we impose a gaussian distribution with mean and variance 1 

// = (0(0)) = 00 , <rl = ((0(0) - 0o) 2 > (3.21) 



as initial condition for P(0, i). We then obtain 

{(j> - 00 + Kt 



21 



2a 2 (t) 



(3.22) 



where 



HH _ gH 



2 £ 3 



Erom Eq. (3.22) it follows that the effect of the linear potential is a mere displacement of the 
center of the distribution. Consequently, the field remains gaussian distributed with mean 
= 4>o — nt and variance <? 2 {t). Nevertheless, to understand the evolution of the distribution 
in terms of the model parameters it is convenient to introduce the timescales for the change 
of the variance a 2 , r CT , and for the change of the mean field //, Tu. These are given by 



where N = Ht is the number of elapsed e-foldings. Also, we have neglected the initial 
variance a 2 to estimate t a . The evolution of P can be described now in terms of r CT and r^. 

When r CT <C the evolution is dominated by the diffusive motion. This implies that the 
coupling between the fields and x is unable to restrain the random motion of the former, 
which then fluctuates as a massless free field. However, since r a increases with time while 
Tfj, remains constant, the diffusive motion of does not persist indefinitely. The quasi- free 
fluctuations of (p can be considered a good approximation to its stochastic dynamics provided 
N <C iVdrift, where 

_ 367T 2 0o H 2 0o r „ 
Adrift = — 3~17 = 17 ■ 3 - 2 5 

Using typical values of the parameters g < 1, e <C 1 and assuming <pQ > H , we conclude that 
the epoch during which fluctuates as a free field encompasses an extremely large phase 
of inflation. Nevertheless, if inflation is sufficiently long-lasting, i.e. for N S> Adrift > the 
diffusive motion becomes subdominant and the distribution drifts towards = 0. Although 
the length of inflation required for that to happen is extremely large indeed, this could be 
the case provided eternal inflation is considered [32]. 

The above result indicates that even in the presence of an effective mass = g 2 (x 2 )iR 
the field manages to fluctuate as a massless, free field. The fundamental feature that allows 
to fluctuate freely is that m 2 ^ is inverse to [c.f. Eq. (3.14)]. The key to understand this 
result physically relies in the heaviness of the x field. As long as m 2 , ~ g 2 4> 2 S> H 2 , the 
response of the system to a random fluctuation increasing (0 2 ) is to adjust = g 2 (x 2 )iR 
to a smaller value according to Eq. (3.14) in much less than a Hubble time. This is to say 
that the swift compliance of the system to adjust (x 2 )ir allows the mean square (0 2 ) to grow 
unimpeded. 



1 The usual initial condition P(<f>,0) = 8{<f> — <f>o) is trivially recovered in the limit o~o = 0. 
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Figure 1. Plot of \xk/Xk\ 2 vs. m x . Solid curves from top to bottom correspond to the time when 
the fc-mode exits the horizon and when the mode exits the coarse-graining scale. The dashed line 
represents the plotted ratio in the limit m x ^> H, i.e. (H/m x ) e 3 . The plot is built using e = 5 x 10~ 2 . 

3.4 Below the crossover scale 

The regime of free fluctuations that we have just described cannot account for the evolution 
of (ft indefinitely. Owing to its random motion, the expectation value of (ft can grow large 
in some coarse-graining patches, thus making \ heavier and keeping our computation valid. 
However, in some patches (ft takes on progressively smaller values. Eventually, the field (ft 
decreases enough so that x ceases to be a heavy field, which then gives rise to the production 
of superhorizon fluctuations of x- I n turn, this production enhances the effective potential 
of (ft making it steeper. The field (ft is consequently driven closer to the origin where x is 
even lighter. As a result of this, the production of superhorizon fluctuations reinforces and 
accumulates as the spectrum of the x field becomes less scale-dependent. Since this implies 
that the <ft is pushed harder towards (ft = 0, the final outcome of the production process is 
that the motion of (ft becomes confined around the origin, thus terminating the phase of free 
fluctuations. 

This trapping of the field (ft occurs similarly to the rapid decay of the inflaton field at 
the end of the broad parametric resonance regime [15], or similarly to the moduli trapping 
mechanism [17] (see also [21, 22]). Nevertheless, in this paper we do not concern ourselves 
with the details of the trapping phase, but rather focus on the ramifications of the existence 
of a scale below which the motion of (ft becomes confined around (ft = 0. In the following we 
refer to this as the crossover scale and denote it by (ft c . We further assume that once the field 
falls below the crossover scale and becomes trapped at the origin, no random fluctuation of 
the system can take the field back to (ft > (ft c . 

From the previous discussion it is clear that the crossover scale is necessarily related to 
the Hubble scale, for x must be a heavy field for (ft to fluctuate quasi-freely. Therefore, it 
seems reasonable to anticipate that setting (ft c large enough so that g(ft c is a few times larger 
than H should suffice to guarantee that (ft performs quasi- free fluctuations when (ft > (ft c - 
However, it is still possible to have a more precise estimate of (ft c . As explained before, the 
crucial feature leading to the diffusive motion of (ft is the inverse dependence (x 2 )ir (ft~ l i f° r 
in that case its effective potential does not become steeper as (ft grows large. This observation 
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constitutes the basis for our estimate, which we explain below. 

We begin by setting fh x = so that m x = g<j). To identify the smallest scale at 
which (x 2 )ir becomes approximately inverse to <p we utilize the scale-dependence of the 
perturbation spectrum of a heavy field (see Eq. (3.15)). Such a scale-dependence implies 
that the quantum average (x 2 )ir on the coarse-graining scale is mainly determined by the 
amplitude of the modes Xk on that scale. This can be clearly appreciated from Eq. (3.14) 
which, apart from constant factors, can be understood as the product of |xfc| 2 ? as given by 
Eq. (3.11), times kg. Using this, we estimate the crossover scale (j) c as the scale at which 
the amplitude \xk\ 2 , rather than (x 2 )ir itself, becomes approximately inverse to m x = gcft. 
In Fig. 1 we depict the ratio \xk/Xk\ 2 , where Xk is the solution to Eq. (3.9) for an exactly 
massless field, for a given k-mode as a function of m x . We plot this ratio at the time of horizon 
exit (upper, solid curve), and when the k-mode exits the coarse-graining scale (lower, solid 
curve). The dashed curve corresponds to the limit of the plotted ratio at very large mass, 
namely 

\Xk/Xk\ 2 ^— e 3 for m x »F, (3.26) 

where we evaluate the ratio at the coarse-graining scale k s . We define the crossover scale 
(f) c as the smallest value of (j) at which the ratio \xk/Xk\ 2 approximately matches (up to the 
first derivative) its value in the limit m x H. Eye inspection of Fig. 1 suffices to estimate 
this scale just by reading m x when the bottom, solid curve and its slope match those of the 
dashed one. This roughly gives 

(/> c ~y/WH/g. (3.27) 

If there is a non-negligible contribution fh x ^ 0, the crossover scale must be set sufficiently 
high so that m x <C g4> c and (x 2 )ir becomes approximately inverse to 4> for eft > 4> c , thus 
allowing the virtually unimpeded random motion of (j). 

We wish to remark that the above is not intended to be an accurate computation, but 
a mere guess based on what is required for the coupled massless field <fi to fluctuate freely. 

4 General expression for field fluctuations 

Having elucidated the growth of fluctuations in the large field regime (<fi 3> H/g), we proceed 
now to describe the probability density in the range 4> > (j> c . As discussed before, if the 
initial distribution peaks around <f> = <f)Q well above the crossover scale, the field manages to 
perform a quasi-free random motion (for N <C N^at)- Owing to this, the expectation value 
of (j) can grow small and approach the crossover scale in some coarse-graining patches. 

As argued in Sec. 3.4, whenever the field's expectation value falls below <j) c its evolution 
becomes confined around eft = indefinitely. It then follows that the probability density 
describing the random motion of <p for (f> > (fr c is not conserved with time, and hence it 
cannot be properly described by Eq. (3.22). To find the appropriate probability density 
function we invoke our assumption on the indefinite trapping of the field around = 0. Such 
an assumption allows us to understand the stochastic dynamics of (j) as that of a brownian 
particle performing a drift-diffusion motion for (f> > (f) c in the presence of an absorbing barrier 
located at cp = <p c . This implies that the probability density is determined by the equation 



d t P = Kd^P+^dlP, 



(4.1) 



with the boundary condition accounting for an absorbing barrier at (f) = (f) c . Such a condition 
is of the form [33-35] 

P(<t> c ,t)=0. (4.2) 

Imposing a gaussian distribution with mean <pQ and variance <Tq as initial condition, the 
solution to Eq. (4.1) for (f> > <p c (where P is positive) is found to be 



P(4>,t) 



l 



\/27r<7 2 
where 

f+ = <Po - 
and C is a constant given by 

C = exp 



exp 



2a 2 



<P- 



8tt 2 k 



C exp 



+ 2<p c — Hit 



(0~^) S 

2a 2 



7T 2 (JqK 



+ 



47T 2 (ToK 

H 3 



(4.3) 
(4.4) 

(4.5) 



The first term in Eq. (4.3) corresponds to the probability density in Eq. (3.22), whereas the 
second term accounts for the loss of probability due to the absorbing barrier at <f> = <j) c . 

The behavior of g~ 1 P((f), t) is illustrated in Fig. 2, where we plot contours enclosing the 
region <7 _1 P(<^>, t) = 10 
(j) c = 3H/g and oq = 

k = W~' i gH 2 in the righthand one to appreciate the effect of the potential tilt k on the 
probability density. The plot shows that for larger values of k, implying a steeper linear 
potential, the field distribution is "pushed" towards the absorbing barrier at an earlier time. 



7 for 7 < 1, 2, 3, 4 (from heavier to lighter shading) using (j)Q = lOH/g, 
5 x 10~ 2 H/g. In the left-hand panel we take k = 10~ 4 gH 2 and 



z 

00 




1.0 1.5 

log 10 (g0/H) 




0.5 1.0 1.5 

logi (g<f>/H) 



2.0 



Figure 2. Contours of g~ 1 P(<j),t) (see Eq. (4.3)). The shaded areas correspond to g~ 1 P(<f> 7 t) = 10~ 7 
for 7 < 1, 2, 3, 4 (from heavier to lighter shading). To illustrate the behavior of P{4>, t) we have chosen 
4>o = 10H/g, 4> c = 3H I g and Co = 5 x 10~ 2 H/g. To demonstrate the effect of the parameter k wc 
choose k = 10~ 4 gH 2 in the left-hand panel and k = 10~ 3 gH 2 in the righthand one. 
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4.1 Fraction of survivors 

In order to trace the approach to the equilibrium distribution we consider the fraction of 
the ensemble in which the field is above the crossover scale, i.e. fluctuating with an out-of- 
equilibrium amplitude, after N e-foldings of inflation. This fraction, which we denote by J 7 , 
can be readily computed by integrating Eq. (4.3): 



P(d>,t) 



l 



[h{<p + ,t)-Ch{y-,t)] 



where 



exp 



2a 2 



1 + Erf 



(4.6) 



(4.7) 



Of course, the fraction of the ensemble in which the field fluctuates with the equilibrium 
amplitude after N e-foldings is l—J-(t). If inflation lasts for a number of e-foldings N <C Adrift 
we can neglect the drift term and consider k — 0. In that case, Eq. (4.6) simplifies to 



J-(t, k = 0)= Erf 



(4.8) 



In Fig. 3 we represent several regions in which T varies within certain range (as indicated). 
Similarly to Fig. 2, we use k = W~ 4 gH 2 in the left-hand panel and k = 10~ 3 gH 2 in the 
righthand one. Again, the plot shows that as the tilt k increases the field distribution is 
pushed towards the equilibrium state at an earlier time. 




o.oi < t < o.i 



o.i < T < 0.95 



r > 0.95 



1.0 -0.5 0.0 0.5 1.0 

logio L?»o-0,)/H] 



1.5 



2.0 




-0.5 



0.0 0.5 1.0 

lOglo [«(0O-0r)/H] 



1.5 



2.0 



Figure 3. Contours of the fraction T varying within certain range (as indicated in the plot). We 
use Co = 0.5H/g and k = 10~ 4 gH 2 for the left-hand panel and k = 10~ 3 gH 2 in the righthand panel, 
where we plot the same contours, to demonstrate the effect of the tilt n. 



4.2 Out-of-equilibrium fluctuations 

Next, we compute the variance ((j) 2 ), which we write as 



/eq 



+ 



(4.9) 
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where ((/> 2 ) e q is the contribution from equilibrium fluctuations, which the field is assumed to 
attain after falling below the crossover scale, and (</> 2 )a is the contribution from the fraction 
of the ensemble in which (p undergoes quasi- free fluctuations. Using Eq. (4.3), the latter can 
be readily computed as 



(j) 2 P(^),t) d(j) 



1 



V27TCT 2 



[I 2 (<p + )-CI 2 (<p-)] 



(4.10) 



where 



poo 


\ (0-^) 2 l 


/ 4> 2 exp 


2a 2 



exp 



2a 2 



o 2 {ip + <j) c ) + 



7T0"' 



(* 2 + 



1 + Erf 



■ (4-11) 



To work with a simple, manageable expression we consider an inflationary phase of 
length N <C ^drift) i n which case we can take k = 0. As previously noticed, this is the 
most plausible case unless eternal inflation, or some other long-lasting inflationary scenario, 
is considered. Using Eq. (4.4) with k = and Eq. (4.10) we find 



exp 



2a 2 J 



a5 + 2(f) c 5 + (a 2 + 5 2 + 



I) Erf 



where we introduce the initial departure from the crossover point as 5 = 
nience. 

If <j) begins above the crossover scale in the entire ensemble (^o S> 



(4.12) 



for conve- 



and al <. 5 2 ) 



then (4> )a should grow linearly with time. Our Eq. (4.12) must then reproduce the epoch 
of linear growth with time for as long as a 2 5 2 . Using that Erf(x 3> 1) ~ 1 and that 
(^ 2 (0)) = a 2 + 4> 2 we find 



! (i))fl * Wc5 + {a 2 + 5 2 + tf) = <</> 2 (0)> + 



4^" 



(4.13) 



which is the expected result [1-3]. Alternatively, this result can be derived from Eq. (4.12) 
simply by taking the absorbing barrier to infinity, i.e. cj) c — > — oo, or 5 — > oo equivalently. 
Next, we obtain the behavior of (<?!> 2 )fi at later times, when the field becomes trapped below 
<f) c in a substantial fraction of the ensemble. Such a situation corresponds to a 2 3> 5 2 . Since 
g(po > H and gcf) c ~ H (see Eq. (3.27)), the condition a 2 3> S 2 implies both a 2 3> 4> 2 and 
a 2 3> cr 2 . Further expanding Erf(|x| <C 1) — 2ir~ 1 / 2 x in Eq. (4.12) we obtain 



l (t)) a ~ 25 



HH 
2^3 



(4.14) 



Comparing this expression with Eq. (4.13) we conclude that once <f> falls below <p c in a 
substantial part of the ensemble the growth of inhomogeneities proceeds at a rate slower than 
in the uncoupled case, with the scaling (</> 2 )fl oc yt. This slowdown is a direct consequence 
of the coupling between <f> and x> which prevents the growth of the field fluctuations (away 
from the equilibrium value) below the crossover scale. 

The behavior of (<^ 2 )fi for several values of the tilt parameter k is displayed in Fig. 4. 
In the left-hand panel we show the case when the field begins slightly above the crossover 
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l°gio £ 2 N log,,, g 2 N 

Figure 4. Plot of (<p 2 ) as a function of g 2 N. We choose C = 3H/g, <r 2 = Q.25H 2 /g 2 and <j> = 3.5H/g 
for the left-hand panel and (f>o = 5H for the righthand panel. The thin curve corresponds to k = 
(Eq. (4.12)) in which the evolution of the distribution is dominated by the diffusive motion. The thick 
curves correspond to the prediction in Eq. (4.10) for three different values of R = n/(gH 2 ). 



scale. We take 4> c = 3H/g, (f>Q = 3.5H/g and a 2 . = 0.25H 2 /g 2 . For these values, a substantial 
fraction of the initial distribution begins with <fi < (f> c . As a result, the probability of finding a 
freely fluctuating field in a member of the ensemble reduces, and so does {4> 2 )& as a result. The 
thin curve represents the prediction in Eq. (4.12) corresponding to a negligible drift motion, 
and constitutes a good approximation to (</> 2 )h as long as N <C iV^rift- For N Adrift 
the motion of the distribution is dominated by the drift, which then suppresses the out-of- 
equilibrium fluctuations. The thick curves in Fig. 4 correspond to the prediction in Eq. (4.10) 
for different values of n (as indicated in the plot). If inflation persists for long enough, the 
fraction of the ensemble with out-of-equilibrium field fluctuations becomes exponentially 
suppressed [c.f. Fig. 3], hence the field fluctuations approach their equilibrium amplitude. In 
the righthand panel we change the initial value to (j>o = 5H/ g while keeping the same values 
for (j) c and ctq. In this case, the field begins above the crossover scale in most of the ensemble, 
which prevents the initial decrease in (^ 2 )fl. 

Note that (</> 2 )fi does not become exponentially suppressed right after departing from 
the behavior with k = 0. Of course, this is because the previous stage of diffusion makes the 
field's expectation value grow large in some coarse-graining patches. Owing to this, forcing 
the field below the crossover scale in those patches is a lengthy process indeed. In turn, this 
implies that the field is able to maintain out-of-equilibrium fluctuations for a longer time in 
the coarse-graining patches where it has grown larger. As a result of this, there may be a 
substantial fraction of the ensemble in which the field has a relatively large expectation value 
compared to the typical equilibrium value. 

To better appreciate and quantify this observation we can combine Eqs. (4.6) and (4.10) 
(see Figs. 3 and 4). Consider, for example, the curve k = 10~ 4 gH 2 in the left-hand panel of 
Fig. 4, which corresponds to (j>o — <p c = 0.5H/g. Using Eq. (4.6) we compute that the field 
is still fluctuating in 1 to 10% of the ensemble when 6 x 10 2 < g 2 N < 4.5 x 10 4 (see Fig. 3). 
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Introducing this lapse of time in Eq. (4.10) we obtain a value of ((ft 2 ) a larger than (ft 2 by a 
factor of order one only. This implies that (<^> 2 )fi is never far away from the equilibrium value 
in a substantial part of the ensemble. 

Also, note that the maximum value g 2 ((ft 2 )^ ~ 38H 2 is attained at g 2 N ~ 5.1 x 10 5 , 
by which time the field is performing out-of-equilibrium fluctuations in 0.2% of the ensemble 
only. Such a result clearly stems from having set the initial field distribution very close to the 
crossover scale, which can be avoided by increasing (fto — <ft c . We consider such an instance in 
the righthand panel, for which (fto — <ftc = 2H/g. From Eq. (4.6) we find that the field performs 
out-of-equilibrium fluctuations in 1 to 10% of the ensemble when 8.7 x 10 3 < g 2 N < 3.6 x 10 5 
(see Fig. 3). Using now Eq. (4.10) in this interval we find a maximum g 2 ((ft 2 )a — 152H 2 , 
which is one order of magnitude larger than (ft 2 (see Fig. 4) and hence substantially larger 
than the amplitude of the equilibrium fluctuations. 

5 Conclusions 

In this paper we have investigated the stochastic dynamics of coupled scalar fields during 
inflation. We assume a massless scalar field 3>, divide it into short wavelength and long 
wavelength parts, integrate out the short wavelength modes and obtain stochastic equations 
of motion for the long wavelength part (ft that we coupled to another (massive) scalar x 
through an interaction of the form g 2 (ft 2 x 2 - We also assume that the coarse grained patches 
thus obtained can be treated as separate universes. The problem is generic but has recently 
been discussed in particular in the context of the flat directions of MSSM [6, 12], and (ft could 
be identified with an MSSM flat direction while % would be a non-flat direction. 

We allow for the possibility that at the onset of inflation, the massless field (ft has a 
large expectation value. As a consequence, the scalar x becomes a heavy field. This allows 
us to integrate out the x field, thus obtaining an effective dynamics for (ft. We find that, 
despite its coupling, the massless field (ft manages to fluctuate as if free. This happens 
because the x field, while being heavy, adjusts the ^-dependent magnitude of its fluctuations 
in a timescale very short compared to the Hubble scale. Owing to this swift response of 
(x 2 ), the random motion of (ft proceeds virtually unimpeded. Thanks to this randomness, 
the expectation value of (ft grows large in some coarse-graining patches and small in others. 
Whenever the field's expectation value falls below the crossover scale, the x field ceases 
being heavy and, consequently, the production of its superhorizon fluctuations is no longer 
suppressed. In turn, the superhorizon fluctuations of the x field enhance the interaction 
potential for (ft. We assume that this enhancement traps the field around (ft = indefinitely, 
while its fluctuations attain their equilibrium amplitude. This assumption is translated into 
the field's stochastic dynamics by imposing an absorbing barrier boundary condition for the 
Fokker-Planck equation located at the relevant scale. 

After integrating out the x field we solve exactly the effective Fokker-Planck equation 
for (ft, thus obtaining an analytical expression for P((ft,t). This is done in subsection 3.3 and 
at the beginning of section 4. Using the result, we compute the fraction of the ensemble in 
which the field's expectation value remains above the crossover scale and the corresponding 
contribution ((ft 2 ){\ to the mean-square. Our expression for (</> 2 )fi reproduces the known result 
((ft 2 ) oc t as long as (ft remains above the crossover scale in most of the ensemble. 

However, when (ft falls below the crossover scale in a substantial part of it, the growth 
of field fluctuations greatly differs from the known result. Instead of growing as ((ft 2 ) oc t, we 
find that ((/> 2 ) fl oc y/t (see Eq. (4.14)). 
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The duration of this epoch of quasi-free fluctuations is nevertheless limited. Beyond 
certain timescale, the diffusive motion becomes overshadowed by the overall drift towards 
<j> = originated by the tilted potential. If the inflation persists beyond such timescale, 
the probability density for field values above the crossover scale goes to zero exponentially 
faster than for a purely diffusive motion (see Eq. (3.22)). Clearly, this result is due, on one 
hand, to our assumption on the existence of an absorbing barrier, and on the other hand, to 
the residual tilt in the effective potential of (j). Therefore, the ultimate fate of the coupled 
massless field <fi is to become trapped around <j> = (where the production superhorizon 
fluctuations of x is n °t suppressed) and fluctuating with the equilibrium amplitude. 

As a final comment, we emphasize that although the number of e- foldings required to 
drive the (f> to its equilibrium state is very large indeed, eternal inflation can obviously result in 
an arbitrary large number of inflationary e-folds for our observable universe. If this turns out 
to be the case, the coupled massless field <p fluctuates with the equilibrium amplitude when 
the scales of cosmological interest exit the horizon in the last phase of inflation. However, 
if inflation does not last for an exponentially large number of e-foldings, the assumed initial 
expectation value of 4> and the subsequent phase of quasi-free fluctuations can give rise to a 
pattern of fluctuations in cf) which can be substantially different from the gaussian pattern 
expected to emerge from equilibrium fluctuations. We illustrate this feature with a particular 
case in which the amplitude of the field fluctuations (^ 2 )h can become up to ten times larger 
than the equilibrium amplitude in 1 to 10% of the ensemble. 
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